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Abstract: Numerous experiments currently underway offer the potential to indirectly probe new 
charged particles with masses at the weak scale. For example, the tentative excess in h — > 77 decays 
and the tentative gamma-ray line in Fermi-LAT data have recently attracted attention as possible 
one-loop signatures of new charged particles. We explore the interplay between such signals, dark 
matter direct detection through Higgs exchange, and measurements of the electron EDM, by studying 
the size of these effects in several models. We compute one-loop effects to explore the relationship 
among couplings probed by different experiments. In particular, models in which dark matter and 
the Higgs both interact with charged particles at a detectable level typically induce, at loop level, 
couplings between dark matter and the Higgs that are around the level of current direct detection 
sensitivity. Intriguingly, one-loop h — > 77 and DM DM — > 77, two-loop EDMs, and loop-induced 
direct detection rates are all coming within range of existing experiments for approximately the same 
range of charged particle masses, offering the prospect of an exciting coincidence of signals at collider, 
astrophysical, underground and atomic physics measurements. 
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1 Introduction 

It has been expected for a long while that new physics would proliferate at about the electroweak 
scale. While the discovery of a Higgs-like particle at the LHC [1, 2] is a profound advance in particle 
physics, there is so far no evidence for other new particles beyond the Standard Model (SM). Although 
new colored particles are strongly constrained, color singlet new particles could still be at large, as the 
LHC has just begun to be sensitive to electroweak processes. It is interesting to explore all possible 
experimental consequences of new light colorless weakly-interacting charged particles. Specifically, 
depending on their couplings, they could possibly lead to the following observations, summarized in 
Fig. 1: 

• Features in the photon spectrum from the galactic center or other astrophysical sources if the 
charged matter couples to the dark matter (DM). 

• Modification of Higgs decay, in particular, h — > 77, if they couple to the Higgs. 

• Electron or neutron EDM if there is CP violation in their couplings. 
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• Modifications of electroweak precision observables if they are weakly charged. Certain parame- 
ters such as oblique parameters are stringently constrained but current bounds on triple gauge 
couplings (TGCs) are still weak. 



Light charged particles could also be produced through electroweak processes at the LHC [3, 4], leading 
to different signals depending on their decay modes. It has been shown by some rough estimates that 
such particles could be within reach of discovery in almost all cases in the 8 TeV run at the LHC, 
and in even the most difficult cases at 14 TeV [3]. However, in this paper, we will not discuss the 
direct detection of these charged particles at colliders. The purpose of this paper is to explore the 
constraints on different indirect observables listed above and possible correlations between them in 
scenarios with light charged matter. In particular, we are interested in two correlations in the cases: 
a) if the charged matter is coupled to both DM and Higgs; b) if there is an order one phase in the 
charged matter sector that cannot be rotated away. 




Figure 1. Possible experimental signals of weak-scale colorless charged matter. 

Now we would like to review briefly all the possible observations and their current experimental 
status. If the charged matter couples to DM, it could mediate loop-level annihilations of DM particles 
into two photons (and possibly photon + Z). This possibility becomes interesting as recently an 
observation has been made of gamma-ray line emissions in the galactic center using 3.7 years of Fcrmi- 
LAT data [5, 6]. While the feature is consistent with one single line at about 128 GeV, a pair of two 
lines at around 111 GeV and 128 GeV gives a slightly better fit to data [7, 8]. (It is unclear if this 
continues to be true in Fermi's updated Pass 8 data, of which only limited information is publicly 
available at this time. It is known that the higher-energy line shifted to 135 GeV after recalibration, 
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and unclear whether the second line can still be detected. It is our understanding that many details 
may still change as the Pass 8 data is validated and analyzed before its public release, so we will take 
a wait-and-see attitude.) The Harvard group further reported possible double line emissions at about 
111 and 128 GeV from unassociated Fermi-LAT sources [9]. These lines could be consistent with DM 
particles annihilating at one loop to 77 (and 7Z for two lines) with considerable cross sections of order 
(crv) ~ 1CP 27 cm 3 s _1 [6, 10]. On the other hand, a similar signal in "Earth limb photons" is difficult 
to reconcile with a dark matter interpretation. In any case, the current cross section at which gamma 
ray lines are potentially detectable is 10 -27 cm 3 s _1 , and whether the tentative lines persist or not, it is 
interesting to consider interpretations of hypothetical new physics at the current sensitivity threshold. 

It is non-trivial to have a DM model give rise to these lines without being constrained by other 
observations. For example, MSSM neutralino DMs annihilating into two photons are ruled out by 
the continuous gamma ray spectrum as at tree-level level, they annihilate into other SM particles 
such as WW and ZZ, whose subsequent decays could also produce photons with continuous energy. 
More specifically, fixing the loop level cross section to diphotons to fit the signal, cross sections of any 
processes contributing to the photon continuum spectrum are constrained to be no more than 5 to 
10 times the loop level cross section [11-13]. Given that, the simplest possibility is that DM is a SM 
gauge singlet annihilating through a loop of light SM charged particles. To get the desired rate, this 
requires quite large couplings and numerical enhancements from coincidences in the mass of DM and 
the particle running in the loop [14-22]. Another possible topology for a rate enhancement is through 
an s-channel exchange of a pseudo-scalar or vector [15, 23-30] with mass tuned close to twice the DM 
mass. The pseudo-scalar (or vector) couples through a loop of light electrically charged particles to 
two photons (or a photon and Z). 1 

If the colorless charged matter couples to the Higgs, it will modify the hyy coupling at the one-loop 
level. Intriguingly, both CMS and ATLAS observed an enhancement in the Higgs decaying to diphoton 
rate while the Higgs decaying to diboson rates a x Br(h — > WW* , ZZ*) are roughly consistent with the 
SM Higgs expectation [1, 2]. If these observations are confirmed by more data, the only explanation 
is to have new charged matter with mass close to the LEP bound, 100 GeV [3, 38, 39]. It should 
be emphasized that these new charged particles should not be colored; otherwise they will be ruled 
out by vacuum instability constraints [3, 40]. For further recent studies on possible mechanisms for 
a diphoton enhancement without enhancing diboson rate, see also Refs. [41-59]. A good, up-to-date 
review of the status of Higgs physics may be found in Ref . [60] . 

On the other hand, if the diphoton enhancement disappears, the data would constrain the coupling 
of a light charged particle to the Higgs, or more accurately, the dependence of its mass on electroweak 
symmetry breaking. 

Furthermore, generically light charged matter coupling to the Higgs can have an order one CP 
violating phase [61], leading at two loops to a signal already constrained by electron electric dipole 
moment (EDM) experiments [62]. In the near future, there will be an update of the electron EDM 
experiment which could potentially improve the current bound by one order of magnitude [63] . Thus 
EDMs could be an independent restriction on possible deviations of the /177 coupling. 

'Two notably different options exist in which the gamma ray lines are actually different shapes with narrow widths 
smaller than the Fermi-LAT current resolution. One is that the lines are narrow box-shaped features from a process 
like DM + DM — > 7Tq + tTq in which 7Tq subsequently decay to 77 and 7Z [31, 32]. The narrow width of the box could 
be explained if 7Tq is degenerate with DM in mass. A simple elegant explanation could be that the 7Tq are in the same 
multiplct with DM due to a symmetry [33]. A UV completion of this scenario was proposed in Ref. [34], The other 
option is internal bremsstrahlung, in which DM annihilating with a t-channel particle emits photons with energy below 
the DM mass [35]. When the mass of the t-channel particle is tuned to be close to the DM mass, the edge could be 
peaked around DM mass [36, 37]. 



- 3 - 



Depending on the quantum numbers, the light charged matter could contribute to both the oblique 
parameters in the electroweak precision analysis and non-oblique parameters such as triple gauge 
couplings (TGC). For the constraints on the oblique parameters, it has been shown in [3, 39] that they 
could be evaded easily. In this paper, we will elaborate more on the constraints on TGCs. 

One topic we will not discuss is muon g — 2. Although formally it fits very well with our theme, 
since magnetic dipole moments are essentially the CP-conserving partners of electric dipole moments, 
phenomenologically its status is somewhat different. The difficulty is that the measured and theoretical 
values disagree, Sa^ « (2.8 ± 0.8) x 10~ 9 [64-66], and this discrepancy is too large to correspond to a 
two-loop effect analogous to the EDMs we discuss, which translate to <5a M on the order of 10 -12 . Hence, 
any two-loop effect is either masked by other new physics that explains the discrepancy (a very exciting 
possibility), or, more plausibly, hidden beneath theoretical and experimental uncertainties that remain 
to be resolved. A one-loop new physics contribution could explain the data, but this requires physics 
like sleptons that carry lepton flavor quantum numbers; for instance, a recent discussion relating g — 2 
to h — > 77 via sleptons appeared in Ref. [67]. We will have nothing to add to its discussion, because 
throughout this paper, we discuss only new charged particles that carry no flavor quantum numbers. 

In Section 2, we will discuss two models that introduce new charged particles running in loops 
to fit a large h — > 77 excess and dark matter annihilation to a gamma-ray line. In both cases, we 
show that loops can induce couplings of the Higgs to dark matter that are in tension with direct 
detection limits. (A third such model is discussed in Appendix A, to avoid tedious repetition in the 
main text.) In Section 3, we discuss two-loop EDM effects and the expected correlation between 
new physics in CP-even and CP-odd observables if phases are generic. (A minor technical detail 
is discussed in Appendix B.) Sections 2 and 3 have some overlap in their content, particularly in 
discussions of corrections to h — > 77 from new vectorlike fermions, but are written so that they can 
be read independently. Concluding remarks arc offered in Section 4. 

2 Induced Higgs-DM coupling via charged matter 

Because the Fermi-LAT gamma ray line and the h — ¥ 77 rate could both be taken as hints of new 
charged matter at the weak scale, it is tempting to postulate new particles that explain both possible 
signals. Indeed, such a suggestion has been made in Refs. [14, 20, 29]. However, when considering the 
full effective theory of the Higgs, dark matter, and new charged particles, one must be careful: couplings 
between dark matter and the Higgs are constrained by direct detection experiments (DM+q — > DM+q) 
as well as (depending on the CP properties of the initial state) indirect detection DM + DM — > h* — > 
WW, ZZ. Even if we initially assume that dark matter and the Higgs are not directly coupled, 
renormalization group evolution in the effective theory will inevitably generate couplings between 
them. The relationship among these various processes is illustrated in Fig. 2. 

2.1 All-scalar effective theory 

First, we will illustrate the idea in the context of a simple model in which the dark matter and the 
new charged particles are all scalars. This is essentially the model studied in Ref. [14, 20], although 
some of our conclusions differ. (The DM/Higgs interactions were discussed much earlier in Ref. [68].) 
We take DM to be a real SM singlet scalar <f>, with = 130 GeV and a Z/2 symmetry tfi — > — <f>. 
The charged matter is a scalar S charged under the SM gauge group as (1, N) Y with N s species and 
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Figure 2. Upper: Possible diagrams contributing to an induced coupling between DM and Higgs. Lower: 
Processes leading to photon continuum and scattering with nucleons at direct detection. The blobs represent 
(one-loop) induced coupling between DM and Higgs through charged matter. 

a common mass ms- The relevant interactions are 

-CD A 0S 2 |5| 2 + X HS \S\ 2 \H\ 2 + \tHcf\H\ 2 + m% a \S\ 2 + X S \S\ 4 + X - m \<\? + X^ 4 

/4|ff| 2 + X H \H\ 4 (2.1) 

with the annihilation processes into two photons depicted in Fig. 3. In the case of more than one 
S species, IS^ 2 should be interpreted as J2i \^i\ 2 an d |£1 4 as (J2i \^i\ 2 ) > one could consider more 
general contractions of the S flavor indices, but there would be no qualitatively different physics. Here 
/j,H = is fixed by the measured Higgs mass rrih ~ 125 GeV, which together with the measured 

Higgs VEV also determines A# s» 0.13. The physical S mass is given by 

m 2 s = m 2 s . + -X HS v 2 . (2.2) 

(In the case that S carries SU{1)w quantum numbers, additional couplings may be present, e.g. 
(W S)(S 1f H) where SU(2)\y indices are contracted within the parentheses. We will not discuss the 
full parameter space of such couplings, which we expect would not qualitatively change any of our 
conclusions.) Thanks to the Higgs low-energy theorem [69, 70], we see that we require Xhs < if 
loops of the S field are to increase the h — > 77 rate. In order to prevent the potential from being 
unbounded from below due to this negative quartic, we require 

A 2 

As > Xs-min = TT ', (2-3) 

allowing a metastable, rather than absolutely stable, vacuum ameliorates this constraint by about a 
factor of 2, according to a tree- level calculation of the bounce action for vacuum decay [40]. 
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DM(0) x ^7 DM(0) % ^ , ' ^jj^' 




DM(0)-'' ^7 DM(0>'' 

Figure 3. The annihilation processes leading to 2 photons. 
2.1.1 The constraint from the gamma ray continuum 

The interaction X^h^IH] 2 provides a dark matter annihilation channel DM + DM — » h* — > VFW, 
with cross section given by [68] 



<™>= X^fcLh/ 1- #7 ^M 2 + — ^? (2 ' 4) 

-«rm* V "V (4m 2 - mU \ m i ) 



i=W,Z 
Xrf,H 

0.028 



3 x lO^cm^s -1 , (2.5) 



taking = 130 GeV and mu — 125 GeV. In the first line = 1 for bosons and 1/2 for Z 
bosons is the Bose factor in the case of identical final state particles. There is also a phase-space 
suppressed annihilation to the hh final state. Note that X^h in our notation corresponds to what was 
denoted A^x/2 in Ref. [14]. Based on studies of continuum gamma rays from the galactic center in 
Refs. [11-13, 71] (as well as comparably strong constraints from radio in Ref. [72]), it appears safe 
to say that an annihilation rate of 10 -25 cm 3 s~ 1 to WW and ZZ is ruled out even with conservative 
assumptions about astrophysical backgrounds, while a slightly more aggressive approach to the data 
would extend the limit down to around 1 to 2 x 10 -26 cm 3 s _1 . We will quote the bound as: 

|A^|<0.05. (2.6) 

Note that different models could shut off this indirect detection channel; for example, Majorana fermion 
DM is in a CP-odd initial state when annihilating, so annihilation through an off-shell CP-even Higgs 
is suppressed. 

2.1.2 Direct detection constraint 

The cross section of the scalar DM <f> scattering off a nucleon through Higgs exchange is 

\^H\ 2 mjf 2 

&SI = i 2 

A ^^5xl0- 45 cm 2 , (2.8) 



0.05, 

where we take the nucleon mass m n = 0.94 GeV and / parametrizes the nucleon matrix element 

(n\m q qq\n) = f g m n [nn], / = ^ = 9 + q E ^ ( 2 ' 9 ) 

q—u,d,s,c,b,t q—u,d,s 
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We use / = 0.30 which is the central value obtained from the analysis in [73]. The most recent update 
of Xenon 100 constrains a si to be smaller than 3 x 10 _45 cm 2 for DM mass at 130 GeV [74]. A similar 
constraint can be obtained for fermionic DM scattering off the nucleon through Higgs exchange. 

In summary, both photon continuum and direct detection constrain the induced DM Higgs coupling 
to be smaller (in absolute value) than about 0.05. 

2.1.3 Matching the data 

The cross section of the DM annihilation to diphotons is 

(H(#^77) = \<***sN 8 r^Aoir^ , (2.10) 

where Qs sums the charge squared over all components inside S and 

Ao(t) = -t + t 2 /(t _1 ) with = m|/m|, f(x) = arcsin 2 ^. (2.11) 

(For Trig < m<p, it is necessary to analytically continue f(x); see, for instance, the discussion of h — > 77 
in Ref. [75].) Demanding (av)(<p<j> — > 77) = 10~ 27 cm 3 s _1 , one obtains X^s^s Qs as a function of 
5 m = ms — TO0. The result is presented in the left panel of Fig. 4. 




Figure 4. Left: Contours of {a{4>4> — > 77)11) = 1,2 x 10 27 cm 3 s 1 in the (ms, X^sNs Yl Ql) plane. Right: 
Contours of ^ 77 in the (ms, XhsNs Qs) plane. 

At one-loop level, DM could also annihilate into jZ, ZZ final states with the first final state 
leading to a second line at m^ — m 2 z / (Am^) = 114 GeV. If the charged matter transforms non-trivially 
under SU(2)y/ : DM also annihilates into WW, which together with the ZZ final state, contributes 
to the photon continuum. Too large a continuum rate relative to photons would be excluded by data. 
It is straightforward to calculate results for general quantum numbers, but we will simply quote the 
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case where the charged scalar S has only hypercharge. The ratio of other annihilation signals to the 
77 line is: 



= 2tan^M/ 1-— 4- «0.4 (2.12) 



ct(77) 

a ( ZZ ) = 2tan 4 g w I 1 - I « 0.006. (2.13) 
cr(77) ^ m; 1 

These formulas can have significant corrections from loop functions if the charged scalar is lighter than 
130 GeV. 

Finally, the charged scalars contribute to the h — > 77 rate 




a x Br(h — > 77) 



M77 " SM 



1 A gs A s £Q 2 ^ A)(r B ) 
2 m| 6.49 



(2.14) 



where t s = Am s /m 2 and - 6.49 is the SM /177 amplitude. We plot Xhs^S Ql as a function of mj 
for /u 77 = 1.5 in the right panel of Fig. 4. 

There are two things to note about Fig. 4. The first is that achieving a reasonable fit both to 
the Fermi-LAT gamma ray line and to an enhanced h — > 77 rate requires both A^s and Xhs to be 
order-one numbers. The second is that, to fit the gamma ray line, it is necessary that ms not be 
much larger than the dark matter mass; otherwise, the coupling needed to achieve a large enough 
cross section rapidly becomes nonperturbatively large. (This raises the intriguing possibility that the 
annihilation process <jxj) — > SS, forbidden today if ms > m<f,, was active in the early universe and 
played a key role in determining the dark matter relic abundance [25, 76].) If we fix a small splitting, 
say 7715 — m = 1 GeV, and consider S to be a set of N$ degenerate states of charge 1, then the 
coupling we need for av = 10~ 27 cm 3 s _1 is already \<psNg — 4.3. (Furthermore, avoiding a potential 
that is unbounded from below requires another large coupling, As > 9.3/iVj.) A 50% enhancement 
in the h — > 77 signal requires A/jsiVg = —2.2. We will now investigate some of the consequences of 
these rather large couplings. 

2.1.4 RGEs 

The one loop RGEs, keeping the scalar quartic couplings, the top Yukawa, and the larger SM gauge 
coupling effects, are presented below, for the case where S is charged only under hypercharge. (See 
related recent work in [77, 78].) 



16t: 2 I3(X h ) 


= 24A 2 ff 


+ \2X H y 2 t + 2A 2 H + N S \ 2 HS - 6yt + \g A 2 - 9giX H 


(2.15) 


16^ 2 /3(A^) 


= 8^1h 


9 

+ 24A«ji ff A4, + 12X H X^ H + 2N s \ HS \ ( f > s + GX^Vt ~ ^I^h 


(2.16) 


16^ 2 /3(A HS ) 


= *\ 2 HS 


9 

+ (4 + 4N s )X HS \s + 4A 0Jf A 0S + 12X H X HS + 6A 0H y 2 - -g%\f H 


(2.17) 


167T 2 /?(A 0S ) 


= 8 A| S 


+■ (4 + 4Ay A 5 A 0S + AXhsX^h + 24A A 0S 


(2.18) 


16tt 2 /3(A s ) 


= (16 + 


4A S )A| + 2A 2 ffS + 2A 2 S 


(2.19) 


167T 2 /?(A ) 


= 72A 2 


f 2X 2 H + N s Xls 


(2.20) 


167r 2 /3(y t ) 


9 3 

= 2 v\ - 


- Sgjyt - ^glvt 


(2.21) 



- 8 - 



(An easy way to keep track of the numerical factors appearing the 0(A 2 ) terms in beta functions 
of quartic terms is to notice that they must compensate the log /i term in the Coleman- Weinberg 
potential, so the beta functions amount to reading off coefficients in Tr.M 4 .) Note, in particular, that 
we have a simple estimate for a coupling between dark matter and the Higgs induced by a loop of S 
fields as shown in the upper left panel of Fig. 2: 

Xhs^sNs A_ X^Ns XhsNs ± log(A/ ms ) 

8tt 2 m<? 4.3 -2.2 N s 2.0 V ' 

Note that log(l TeV/ms) ~ 2, so the log will already have this size even when running from quite 
a low scale. It is apparent that our bound \\$h\ ^ 0.05 from direct and indirect detection is in 
some tension with our desire to explain both the Fermi-LAT gamma ray line and an enhancement in 
h — > 77. The problem is ameliorated when the number of species, Ns, is large, but this also makes 
the renormalization group effects large. 

For N$ = 1, an even more immediate problem is that avoiding a potential unbounded from below 
requires A5 9. Then the leading term in the As beta function is y^^OA 2 , « 11, and there is no 
sense in which the theory is under perturbative control. For a single charged scalar (of charge 1), 
it is simply not possible to discuss a 50% enhancement in h — > 77 while maintaining a perturbative 
theory. Hence, we should focus attention on Ns > 2. The large couplings suggest that our RGEs 
will become nonperturbative at low scales. We can quantify this by defining, for each coupling, a 
perturbativity limit at which the beta function becomes > 1 when all other couplings are turned off. 
(This definition was used in Ref . [78] . It usually corresponds to smaller couplings than those for which 
the two-loop beta function is larger than the one-loop beta function, but in practice we find that a 
coupling exceeding this bound will very quickly run large enough to exceed the other as well.) For 
example, we define A™ ax by the condition T ^24(A^ ax ) 2 = 1, i.e. A£ ax = ^/2jZir ~ 2.6. For each 
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coupling, we define an analogous A™ ax and define a normalized coupling by A^ = Aj/A™ ax . 

We plot some examples of RGE evolution for the normalized couplings in Figure 8. We begin the 
evolution by fixing A# from the Higgs mass, A^s-^S = 4.3 to fit the gamma-ray line, As = 9.3/iVj 
for vacuum stability, XhsNs = —2.2 for a 50% enhancement of h — > 77, y t from the top mass, and 
Xcf>H = -V — at low scales. The result is that, for Ns = 2, perturbativity is lost almost immediately 
on evolving to higher scales, while for Ns = 3 it is lost around the TeV scale. This suggests that the 
most reasonable interpretation of these models is as composite models, where the scalars are bound 
states. This makes large couplings As and A^s seem more natural; however, the large value of Xhs 
suggests that perhaps the Higgs would be composite too. In this case, the small values of X^h (required 
by direct and indirect detection bounds) and Xh (required by the Higgs mass measurement) seem hard 
to reconcile with the idea that composite states would generically be strongly coupled to each other. 
We could move to larger values of Ns to postpone the loss of perturbativity to higher energy scales, 
but would still face a puzzle in the small value of X^h- For instance, by choosing Ns — 6, we can 
postpone the loss of perturbativity to a scale of 250 TeV. Then, Eq. 2.22 has a factor of log — ; we 
have increased both the log and Ns, and so have not really helped solve the problem. Indeed, solving 
the RGE shows that A^s blows up first; the RGE I671" 2 /3(A^s) = SA 2 ^ tells us that the coupling will 
blow up at about 



A uv 2tt 2 

log « — 

ms X<f,s 



^ (2.23) 
4.3 V ' 



Thus, the factor of jjj- log in Eq. 2.22 does not scale with Ns, and the problem is unavoidable: 
generically, RG evolution is in conflict with the lack of a direct detection signal. In short, although 
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Figure 5. Renormalization group evolution of "normalized couplings," with the number of scalar species 
Ns = 2 at left and 3 at right. The normalized couplings A,; = A^/Af 13,51 , where A™ ax is determined by 
/?(A™ ax ) = 1 even when all other couplings are set to zero. Hence, when any Xi = 1, we expect that perturbation 
theory is no longer reliable. The dashed horizontal line at 1 is to guide the eye to the approximate perturbativity 
boundary. For Ns = 2 this happens almost instantaneously, while for A^2 = 3 it happens at about 1.4 TeV. 

it is very appealing to consider the idea that new charged particles explain both dark matter annihi- 
lation and the h — ¥ 77 enhancement, in this simplest model our closer look has shown that quantum 
corrections spoil the nice idea. 

Aside from raising N$, one could try to ameliorate the Landau problem by making the electric 
charge of S larger. For instance, making S charge 2 instead of 1 will reduce the required values of 
X^s and Xhs by a factor of 4. This brings the loop induced X^h to about the level of the XENON 
bound provided the cutoff scale remains at about 1 TeV. Note that S, and any other new charged 
particles we discuss in this paper, must decay promptly to avoid strong collider bounds on heavy 
stable charged particles (HSCPs). At large charges of S, any effective operator allowing it to decay 
to charged standard model particles will be of high dimension, so that larger charges rapidly lead 
to longer lifetime decays (or require extending the model with a bevy of new intermediate states of 
progressively smaller charges). A detailed analysis of S decays and collider constraints is beyond the 
scope of this paper, but we expect that our conclusions can be weakened only mildly by constructing 
models where S has larger charge and has escaped detection so far. 

Finally, we comment that the analysis of Ref. [20] is rather similar to ours, including also RG effects 
and observing the existence of a low-scale Landau pole. However, one of our central conclusions, that 
attempts to fit the Fermi-LAT line and the h — ¥ 77 rate are necessarily in tension with direct detection 
bounds, is less clear in their analysis. Their discussion of thermal relic abundance, which we have not 
considered here, is interesting. 
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2.2 Resonant annihilation 

A scenario that has received a great deal of attention as a possible explanation of the Fermi-LAT 
gamma line is resonant annihilation [15, 23-30]. By exploiting a pole in a propagator to enhance the 
annihilation cross section, the large couplings and nearby Landau poles we encountered in the previous 
subsection may be avoided. The cost is tuning the mass of an intermediate particle to be close to twice 
the dark matter mass, for no deep reason. We will consider primarily the case of Majorana fcrmion 
dark matter in this subsection. Two identical Majorana fermions, annihilating at low velocities, form 
a CP-odd initial state, so the intermediate particle should be a pseudoscalar rather than a scalar. (In 
other words, an intermediate scalar would lead to a suppressed p-wave, rather than s-wave, annihilation 
cross section.) This might seem to preclude the generation of Higgs-related signals, since the Higgs 
does not mix with pseudoscalars in the absence of CP violation. However, any new charged particles 
that affect the h — > 77 rate would, in general, allow the presence of new CP-violating phases, in the 
presence of which a direct detection signal could be generated. Alternatively, although we will not 
discuss it in detail, one could consider Dirac fermion dark matter, which could annihilate through a 
CP-even scalar rather than a pseudoscalar. 

2.2.1 s-channel exchange of a boson 

We will take DM to be a Majorana fermion £ with mass toj, annihilating through an intermediate 
pseudoscalar boson (f> with mass and width T^. The intermediate state <fi couples to photons 
through AL, species of charged light fermions ip and tp c (with conjugate charges) with mass and 
charge Q^. The Lagrangian is given by 



1 i 



(2.24) 



(with g^ , g^p real couplings and <fi one real degree of freedom) leading to the annihilation process as 
depicted in Fig. 6. We will take the dark matter mass mj = 130 GeV. 



DM(£) 



DM(£) 





1> 



Figure 6. The annihilation processes leading to two photons in resonant annihilation model. 

The formula of the cross section for ££ — > 77, allowing for multiple fcrmionic states to run in the 
loop, is given in [25, 79, 80]: 



47T 3 



( 4m i - m i) 



^2 n *Q%9tI>™">I>I 




(2.25) 



where for m$ > m s , I (^|) = - f arctany ^^^ j . Demanding (aV)(f£ 77) = 10 27 cm 3 s 1 , 

and assuming all the fermions in the loop are degenerate, one obtains g^g^pN^Q^ as a function of 
S m = m^j, — m^. We determine by summing the partial widths <fi — > ££, <fi — > ipip, and — > 77 
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(when the modes are kinematically accessible) using the formulas in Ref. [25]. Because the partial 
widths depend on g^ 7 g^p, and the charge and species count of tp independently, we plot only the case 
5? — 9-4) = = 1. The result is presented in Fig. 7. We see that the tp particles should not be 
much heavier than £ if we want to fit the line without large couplings. More importantly, we require 
an approximate resonance condition, 2m^. With a high degree of fine-tuning the couplings can 

be made quite small; e.g., we can achieve g% = g^ « 0.1 for \m^ — 2m^ \ « 0.5 GeV. With less extreme 
fine-tuning, e.g, a 10% coincidence in masses, the required couplings are approximately 1. 
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Figure 7. Resonant annihilation. Left: contours of (o-(££ — > yj)v) = 10 -27 cm 3 s _1 in the (Sm,g^ = g,p) 
plane, where Sm = m^ — m^. Purple, solid: = 270 GeV; red, dashed: m^, = 300 GeV. At right: the same 
contours as a function of m^; purple: Sm = 50 GeV; red, dashed: Sm — 2 GeV. 

R we also want to alter the h — » 77 rate through loops of this charged particle, we need a more 
complicated structure. Following Ref. [3], we assume that tp carries only hypercharge while another 
field x carries SU(2) W charge, tp, tp c ~ (1, l) T i, X, x c ~ (1> 2)±i ■ The Lagrangian is: 

- C = m^ipip c + m xX x c + yHipx + y c H^tp c X c + cc. (2.26) 

and we also add, for generality, a coupling g x (f>x c X- We will be interested in the general case with 
a non-zero CP violating phase arg(j/j/ c m^m* ) ^ 0. At loop level, this will also force us to consider 
complex values of the couplings g^ and g x . In this case, the one- loop RGEs for our new Yukawa 



- 12 - 



couplings (neglecting SM gauge coupling effects) are extracted from the general formulas [81-84]: 

167r 2 /3(g € ) = % (2|. 9e | 2 + |^| 2 +4|.9 x | 2 ) (2.27) 
16tt 2 /3(^) = H (3 M 2 + |. 9? | 2 + 4 \g x \ 2 + \y\ 2 + |y«f ) - 4<?t yy c (2.28) 
16^ 2 /3( 5x ) = g x (5 \g x \ 2 + \g s \ 2 + 2 \g f \ 2 + \ (ky| 2 + [yf - 2g\ytf (2.29) 
16^ 2 /?(y) = \y (l5x! 2 + + 5 \y\ 2 + 2 |y c | 2 ) - 2y^ 9i ,g x (2.30) 
16^ 2 /3(y c ) = \y c (| 5x | 2 + + 5 \y c \ 2 + 2 |y| 2 ) - 2yt 9lP g x (2.31) 

We will evolve these RGEs to examine the extent to which annihilating through resonant enhancement 
relieves the Landau pole problem of the all-scalar model. 

The modified h — > 77 rate, in the limit of large m^, iX , is determined by the low energy theorems 
to be [3, 61, 62, 70]: 



M77 



i+X-l^flT— logdetA^M 
.Asm 3 ^Slogv 



2 



2 , d 

a — Q % m — arg det M 

Asm ^ dXogv 



■2 



(2.32) 



where the first term arises from the modified CP-even hF^F tJ,l/ vertex and the second term from 
the CP-odd h,F jlu F llv term. Here _4sm = —6.49 represents the SM amplitude. As in the scalar case, 
the result is modified by familiar loop functions that correct for finite mass, which are given in the 
CP-even and CP-odd case respectively by [85] 

A 1/2 (r f ) = ~ T/ + T f ) arcsin 2 ^\ (2.33) 

A1/2 (Tf) = T f arcsin 2 y^, (2.34) 

(with Tf = 4m 2 r/m 2 ) and asymptote to 1 when 2m/ 3> rrih- (Note that there is a minor error in 
Ref. [86], which assumes the same loop function for the scalar and pseudoscalar decay modes.) 

We choose two representative points in parameter space that fit a 50% enhanced h — > 77 rate. In 
both cases, we arrange for the light mass eigenstate to be at 140 GeV, such that it is near the DM 
mass but slightly too heavy for the dark matter to annihilate into two of our new charged fermions. 

• Moderate phase: m^, = m x — 346 GeV, y — 1.37e - 27 ™, y c — 1.37. The mass eigenvalues are 140 
and 577 GeV. As we will see in the next section, this point is excluded by the nonobservation 
of an electron EDM: it predicts d e /e = 9.3 x 10~ 27 cm. However, this exclusion can be avoided 
if the EDM is canceled by another contribution which must be tuned to about the 10% level. 
(And note that our phase is already somewhat small; maximal CP violation would require the 
EDM to be tuned at a few-percent level.) 

• Small phase: — m x — 333 GeV, y — l.lle ' 027 ™, y c = 1.11. The mass eigenvalues are 140 and 
526 GeV. This point is currently safe from EDM constraints, predicting d e /e = 7.3 x 10~ 28 cm, 
but would likely be detected with next-generation electron EDM measurements [63]. 

In both cases, we will also consider equal couplings g = g^ = g^ = g x of the pseudoscalar resonance, 
chosen to achieve a gamma-ray line rate of 1.0 x 10 -27 cm 3 /s. These turn out to be relatively insensitive 
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Figure 8. Renormalization group evolution of "normalized couplings" for the case of annihilation through a 
scalar resonance. Parameters are as in the "moderate phase" case described in the text; the "small phase" 
case has similar RG evolution. At left, we have chosen m^> w 300 GeV, far enough off-resonance that the 
Fermi- LAT line rate requires g = g% = g^ = g x w 1.1. In this case the couplings become non-perturbative at 
scales of order 100s of TeV. At right, we show results for « 270 GeV, closer to the resonance with g ~ 0.55. 
This significantly increases the perturbativity range for the couplings g. Thus, unlike the all-scalar model for 
the gamma ray lines, in this scenario the larger loop contributions arise from fitting the h — > 77 rate, and in 
fact the Higgs quartic running (not shown here) is the most important effect that will require new physics at 
low scales [3]. 



to the mass of the heavy eigenstate. We consider two choices of resonance mass: = 270 GeV, which 
requires g s» 0.55, and = 300 GeV, which requires g rj 1.1. The renormalization group evolution 
of these couplings is plotted in Figure 8. The main qualitative feature to note is that, compared to 
the model in Section 2.1 without resonant enhancement, the couplings stay perturbative until much 
higher scales. There is another crucial RG effect not visible in our plots, discussed in detail in Ref. [3]: 
the large Yukawa couplings y and y c will drive the Higgs quartic coupling negative at low scales. This 
vacuum instability suggests that the model must be altered, for instance by adding light superpartners 
of the ip and \ fields to cancel their effect on the Higgs quartic. However, we should emphasize that 
if one only wishes to explain a gamma ray line signal, this concern is not relevant, and the couplings 
involved in dark matter annihilation can remain perturbative to very high energies if the (f> mass 
is tuned so that annihilation is on resonance. 
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Figure 9. The induced Higgs/pseudoscalar mixing through a CP-violating loop of tp an d X fermions. 
2.2.2 Direct detection 

In the presence of CP violation, scalars and pseudoscalars can generically mix. Integrating out ip and 
X, we find that the Coleman- Weinberg potential contains a term, generated by the loop in Figure 9: 

lm(vv c ) I m?i, m?i, frit mt\ . 

V CW ^ -^L± ^ + g ^ mx) ^ + _±_ log _ log _x j flt^. (2 . 3 5) 

when expanding around the origin. Here we have taken the mass terms my, , m x to be real, so that 
the phases would be contained solely in the Yukawas y, y c . A nontrivial phase can lead to significant 
mixing. To calculate the mixing, we expand the Coleman- Weinberg potential around the physical 
Higgs VEV and evaluate at a renormalization scale fx = ^Jm^m^. We tabulate the resulting mixing 
angles between Higgs and (j), 9h4>, for four cases (moderate or small phase, on or off resonance) in 
Table 1. 

The mixing between <f> and higgs lead to an effective operator Dj^Dqq, with the Dirac fermion 
D = (£,J<t 2 £*) t and q the SM quarks, relevant for direct detection. In the non-relativistic limit, one 
could see that the scattering rate off nucleons through this operator is suppressed by the momentum 
transfer (the amplitude is oc ■ q), leading to an interesting recoil spectrum rising at high recoil 
energy [87]. However, the small rate limits its possibility of being detected. The current XENON 
direct detection limit is g^6h<j> ^ 15(771^/300 GeV) 2 , about two orders of magnitude above g^Ohcf, ~ 0.1 
in the case with a moderate CP phase. Similarly, the conjunction of pseudoscalar and scalar couplings 
leads to a r>wave suppressed indirect detection cross section for annihilation through the Higgs/0 
mixture. Thus, the pseudoscalar resonance models are safely out of reach from direct and indirect 
detection bounds from loop induced DM-Higgs couplings for the foreseeable future. 



dh<jj 


Phase 0.2vr 


Phase 0.02tt 




= 270 GeV 


0.076 


0.0047 


TO0 


= 300 GeV 


0.11 


0.0068 



Table 1. Mixing angles between Higgs and the pseudo-scalar (j>, for both moderate and small phase and on- 
or off-resonance. 
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3 Correlation between CP-odd and CP-even observables 



3.1 Operators and corresponding observables 

We will assume that new charged matter does not interact or mix with the SM fermions at tree level 
and they do not contribute to the EDMs of the SM fermions at one-loop order. (Thus, we will also 
not discuss the discrepancy in the measured muon g — 2, which would typically require new physics 
with leptonic quantum numbers exerting a one-loop effect.) The one- loop EDM is generically ruled 
out unless the new CP violating phases are tuned to be small (< 10~ 2 ). However, new charged 
particles could still contribute at the two-loop order to the SM fermion EDMs through Barr-Zee type 
diagrams [88]. To see the correlations between the CP-odd observables (including EDMs) and CP- 
even observables, it will be useful to perform an operator analysis first, which strictly speaking is only 
valid in the limit when the charged matter is heavy and could be integrated out. We will use operator 
analysis to clarify the correlations of the observables, whereas for the numerical evaluations we will use 
the full-fledged loop calculations. Charged matter with physical phases contributes to 6 CP violating 
dimension-six operators built out of the Higgs and the SM gauge fields. Among them, two involve 
the SU (3) color field strength, one of which is the famous Weinberg operator [89] . We only inspect 
the four operators generated by loops of colorless particles. They and their corresponding CP-even 
operators with similar structures are, following the notations in [90], 

Ow = e ahc W« u W h v x Wf , O w = e abc W^W* x W c x » 

O hW = H^HW^W a ^, O hW = H^HW^W a ^ 

O hB = H^HB^B^, O hB = rfHB^B^ 

O wb = (H^ a H)W^B^, O wb = (H^ a H)w;^, (3.1) 

where a a denotes the three Pauli matrices and a is the isospin index. The operators have coefficients 
dj bounded by the interval [— 1/A^ eg , l/Ap OS ], where A is some high scale. 2 

Now we specify the observables these operators contribute to. ft is well known that Owb gives 
the S parameter in the electroweak precision tests (EWPT), 

_ 4swcwv 2 awB 
a 

where sw = sin#iy, cw = cos 9w, and v = 246 GeV. 

Among all the operators, O w ,0 W b {O w ,O hw ,0 WB ) modify CP-even TGCs (CP-odd TCGs). 
More concretely, the general triple gauge couplings could be parametrized as [93] , 

Cwwv/9wwv = igj (w+ w-<*v v - h.c) + inyW+W-V^ + l**w+w- vX v£ 

m w 

+ik v W+W-V^ + ^ W +W-» X V£ + • • • , (3.2) 
m w 

where V is either Z or 7 and gwwz = — e cot 9w , 9ww-y = — e. The first line of Eq. 3.2 contains 
CP-even TGCs while the second line contains CP-odd TGCs. The dots represent C-violating TGCs 
arising from operators at high orders in the SM effective theory. In the SM, g\ = 1, n,y — 1 and 



2 In some literature [91, 92], more CP-odd operators were listed. As we show in Appendix B, those additional operators 
can be written in terms of the operators in Eq. 3.1 using the equations of motion. 
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Ay = at tree level while kv and Xy are zero even at one-loop order in the SM due to unitarity of the 
CKM matrix. The contributions to the parameters in Eq. 3.2 from high-dimensional operators are 

_ v 2 s w _ v 2 c w 

OKz — OKy — d-WBi 

Cw Sw 

OAz = 0A 7 = , 

9 

c~ o 2 , v2s W ,~ „ 2 "^ff 

6k z = 2v a hw H a WB , 6k 7 = 2v a h ^ a^ B , 

Cw S\y 

8Xz = 8X 7 = 6m w aW . (3.3) 

The operators O hw , O hB , W b, O hw , O hB , O wb also modify the Higgs decays 
r(h -> 77 ) 



S-r-isl ii, „ — r 1 . in ,. _ s „ r„ n„ ,.\ ' 8nV 2 (Sy i r0. h w + <^yd h Q — SwC\V^Wb) |-> |, 



1 , 8™ (s w a h w + c w a hB - s w c w awB) 



aAs 



M 



M 



where .Asm = Aw + -At ~ —6.5 is proportional to the SM amplitude. 

The four CP-odd operators also contribute to the electron EDM at the one-loop order through 
Barr-Zee type diagrams [88], 



Of 
e 




SwCw^WB m I T2 



A 2 

WB , 



(3.5) 



where T3 is the isospin of the charged matter. 

Now we review the experimental status of the observables. Currently the most stringent constraints 
on CP-even TGCs are from the measurement of differential cross sections of e + e~ — > W + W~ at LEP 
2 [94] and those of pp -> W + W~/WZ -> lujj at CMS [95]: these lead to constraints -0.038 < \ Zn < 
0.03, -0.11 < Ak 7 < 0.14 at 95% CL. ATLAS also measures CP-even TGCs in the fully leptonic 
channel, obtaining weaker constraints due to the smaller branching fractions [96]. The CP-odd TGCs 
are also studied at LEP and the Tevatron [97, 98]: k z = -0.09±j^, X z = -0.08±0.07, |A 7 | < 0.32 at 
95%. The Higgs data start to constrain the modifications of Higgs couplings, though due to the large 
uncertainties the constraints are weak at the moment. Nonetheless, it is very interesting that Ohw 
and OhB are becoming as important as the traditional electroweak precision operators. For EDM 
experiments, the current electron EDM bound is d e /e < 1.05 x 10~ 27 cm (90%) [99, 100] and the 
neutron EDM is also stringently constrained: d n /e < 2.9 x 10~ 26 cm (90%) [101]. (A bound on the 
EDM of mercury, d( 199 Hg)/e < 3.1 x 10~ 29 cm (95%) [102], is also noteworthy.) It is expected that 
there will be an update in the electron EDM measurement in the near future, improving the current 
bound by an order of magnitude [63]. We summarize the current experimental constraints on the 
coefficients of these operators in Table 2. 3 



3 One should be careful in interpreting the bounds. For example, O^w^hB are always generated with a coefficient 
OEM /t ~ 10 -3 . Thus Higgs data itself does not put a strong constraint on the mass of the charged particle so far. 
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o 


Experiments 


Ai(TeV) 




EWPT [103] 


12.6 (90%) [104] 


OhW, OhB 


h -> 77 [105, 106] 


4.1 (68%)/ 7.4 (68%) 


Ow 


CP-even TGCs [94, 95] 


1.3 (95%) 


®WB 


CP-odd TGCs [97, 98] /electron EDM [99, 100] 
CP-odd TGCs [97, 98] /electron EDM [99, 100] 

electron EDM [99, 100] 
CP-odd TGCs [97, 98] /electron EDM [99, 100] 


0.5 (95%)/ 38 (90%) 
0.9 (95%)/ 24 (90%) 

48 (90%) 
0.5 (95%)/ 35 (90%) 



Table 2. Current experimental bounds on operator coefficients (the CL are in the parenthesis). The operator 
coefficient a, is bounded by the interval [— 1/A™, l/Ap OS ]. The Aj (in TeV) shown in the table is the average 
of A neg and A pos . When en is experimentally bounded to be negative (positive) definite, we only quote A neg 
(Ap s). For constraints from h — > 77, we fix Higgs mass at 125 GeV. 



3.2 Correlation between CP-odd and CP-even observables 

Now we want to explore possible correlations between CP-even and odd observables, in particular, 
the correlation between Higgs observables and EDM experiments. To enhance the Higgs diphoton 
coupling, electroweak symmetry breaking needs to contribute negatively to the charged matter mass. 
This can be realized without introducing any physical CP phase, for example, through a single scalar 
with a large negative quartic coupling — AslSpl-ff | 2 and A5 > as already discussed in Sec. 2.1. 
More generally, however, this is realized with presence of new CP phases, for example, in models with 
vector-like matter fields which obtain part of their masses from electroweak symmetry breaking. The 
general mass matrix, e.g., for fermions, is 



Cm = -(^ Q X +Q ) I! * \ r_ n +cc, (3.6) 




with the Higgs VEV given by (H) — v/y/2 = 174 GeV and ip,x are Weyl fermions. There is one 
physical phase, <fi = arg ^m^,m* yy c ^j , that cannot be rotated away by field redefinitions. In terms 
of operators, the diagrams generating CP-even operators, Ow , Ohw, @hB lead to O h yy,O h g with 
insertion of the physical phase. Notice that the WWW operator is not generated at one-loop. The 
reason is that the W's and the Z only couple to fermions of the same chirality. Without Higgs 
insertions, as each mass insertion flips chirality, the diagram is always proportional to even powers of 
|m^,| 2 or |m x | 2 , which are always real. The WWW operator could be generated at the two-loop order 
or, similar to the Weinberg operator GGG, WWW receives a finite threshold correction from a heavy 
SU(2)w charged particle with a non-zero EDM d e and mass m, 

where g is the SU(2)w coupling. The constraint on a w from EDM translates into 



T 3 m 



< 5x10 



lToV 
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3.3 EDM 



As shown in Table 2, EDMs are more powerful CP-odd observables compared to the CP-odd TGC 
measurements. The current bound on the electron EDM is d e /e < 1.05 x 10~ 27 cm, which will be 
improved by an order of magnitude in the near future [63]. If a non-zero electron EDM is observed 
then, there are three possibilities: 

1. One- loop EDM if the charged matter has lepton quantum numbers with small CP phases (< 
10- 2 ); 

2. Higher-order contributions that will also affect the Higgs decays: in terms of high-dimensional 
operators, O h] ^,O h ^ and their CP-even counterparts are generated. This could originate from 
vector-like charged matter with chiral mass terms; 

3. Higher-order contributions that do not affect the Higgs decays, e.g., only WWW is generated, 
perhaps from a heavy weakly-charged particle with a non-zero EDM as discussed in the previous 
subsection. 

In this section, we focus on the second possibility. We evaluate the EDM bounds and reach using 
the two-loop EDM formula in [91] and the modification of Higgs decaying to diphotons in two concrete 
models as discussed in [3]. The two models are 

Vector doublets + singlets ("vector-like lepton"): ip,ip c ~ (l,2) ± i, \, \ c ~ (1, l)=pi ■ The 
Lagrangian leading to (3.6) is 

- L = m^i/S + m x xx c + yHtpx + y c H^ c X c + cc. (3.8) 

Vector doublets + triplet ( "wino-higgsino" ) : ip, ip c ~ (l,2) ± i, x ~ (1>3)q. We identify x 
and x°\ the Lagrangian leading to (3.6) is 

- L = m^ijj + ^m x xx + V2yHipx + V2y c H^ c X + cc. (3.9) 

In both cases, we define <j> =arg(yj/ c m^,m* ) € (0, tt/2) and N the number of species for these fermions. 

The results are shown in Figs. 10 and 11. To evade the stringent EDM bounds and have a 
diphoton enhancement /i 77 > 1.5, the CP-violating phase has to be small: <fi < 0.037T for the "vector- 
like lepton" model or cf> < 0.057T for the "wino+Higgsino" model. Besides, for larger </>, it is difficult 
to get a diphoton enhancement /x 77 as the enhancement mainly comes from the interference between 
the real part of the new fermion contribution and the SM contribution (for N = l). 4 For fixed 
light and heavy mass eigenvalues, the real part of the new fermion contribution, which scales with 
\ljy c \ = yy Ccos 4>i decreases as <fi gets bigger. Thus if the diphoton enhancement is confirmed, it will 
pose a new interesting Higgs CP problem at the weak scale analogous to the SUSY CP problem: Why 
do those new charged degrees of freedom that modify Higgs decaying to photons have 
small phases? 

Neutron EDMs lead to comparable constraints. The case where new charged particles are scalars 
(like staus) rather than fermions leads to a somewhat different story. Even if there arc nontrivial 
phases in their mass matrix (e.g., in the A-terms in the case of staus), integrating out such particles 
does not generate the operator hFF. Although this can be straightforwardly seen from the loop 

4 For N > 2, the contribution from the imaginary part itself could be important as it scales as N 2 while the interference 
between the real parts scales as JV. 
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Figure 10. Upper: "vector-like lepton" model; Lower: "wino-Higgsino" model. N = l,m^ = m x ,y — y c in 
all these plots. cj> =arg(j/i/ c m^m* ). The horizontal and vertical axes correspond to the light and heavy mass 
eigenvalues. The solid purple line is the current EDM constraint d e /e — 1.05 x 10~ 27 cm with the grey region 
excluded; the dashed purple line is the projected constraint d e /e = 10 -28 cm. The green lines denote the 
diphoton enhancement ^t 77 . 

calculation, it can also be understood as a consequence of the fact that the argdet A4 coupling arises 
from an anomalous rotation of fermion fields, whereas scalars have no anomalies. However, if there 
is a pseudoscalar particle in the spectrum that can run in the two-loop EDM diagram in place of the 
Higgs, or if CP-violation leads the Higgs to have a small pseudoscalar-like coupling to the electron 
(e.g. by mixing with a pseudoscalar), there will still be a two-loop EDM [107]. Thus, in the case of 
charged scalars, the Higgs CP problem would be less robust: if all pseudoscalars are heavy, the EDMs 
can be rather small. (There are other difficulties for such an interpretation of an increased h — > 77 
rate, as new charged scalars typically have vacuum stability problems [40, 46], although there is still 
viable parameter space for quite light scalars [55].) 
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Figure 11. Upper: "vector-like lepton" model; Lower: "wino-Higgsino" model. N = 2, my, = m x ,y — y c in 
all these plots. The horizontal and vertical axes correspond to the light and heavy mass eigenvalues. The solid 
purple line is the current EDM constraint d e /e — 1.05 x 10~ 27 cm with the grey region excluded; the dashed 
purple line is the projected constraint d e /e = 10 -28 cm. The green lines denote the diphoton enhancement 

4 Conclusions 

Charged weak-scale matter is the key ingredient in most explanations of h — ¥ 77 enhancement and a 
monochromatic photon line at around 130 GeV in the Fermi data. It is tempting to have a unified 
explanation for both of them in which some charged matter couples to both Higgs and DM. We find 
that a fine-tuning at 10% level or worse is inevitable in these models. The large couplings required 
by the h — > 77 enhancement and the photon flux responsible for the line would lead to low-scale 
Landau poles. Even worse, they could induce a considerable DM-Higgs coupling radiatively even if 
this coupling vanishes at tree level. To evade the direct detection and photon continuum constraints, 
one needs to invoke a 10% level tuning between the tree-level and radiative DM-Higgs coupling. In 
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certain models such as the resonant annihilation models, the large couplings and Landau poles could 
be avoided by tuning the mass of the intermediate boson to be close to twice the DM mass for no 
deep reason. For example, with a 10% coincidence in mass, the required couplings are approximately 
1. Constraints on DM-Higgs couplings from direct detection could be avoided in models with DM 
as Majorana fermions and small CP violating phases. In light of these and other difficulties, like the 
vacuum stability problems that plague h — > 77 models, it seems likely that one or both signals will 
fade away; the alternative is not only that light charged particles exist, but that a spectacular array of 
interesting new physics will be found at nearby energies. A more optimistic way of stating this result 
independent of observed anomalies is that the current experimental sensitivity to h — > 77, gamma-ray 
lines from dark matter annihilation, and dark matter couplings to the Higgs induced at one loop by 
charged weak-scale particles are all comparable, so that if new particles are lurking just below the 
current limits, we can hope that a spectacular set of correlated signals will emerge. 

A new CP violating phase could also be present in the couplings of new charged weak-scale 
matter. Even if the new matter does not couple or mix with the SM fermions at tree-level, they 
could generate two-loop EDMs through Barr-Zee type diagrams. In particular, in models with new 
vector- like fermions to enhance the diphoton rate to be more than 1.5 times the SM rate, an order-one 
new CP phase leads to too large an EDM that is ruled out already unless a cancelation with other 
contributions, corresponding to a few percent tuning, is present. This will impose an interesting Higgs 
CP problem, that is, if the diphoton enhancement is true and caused by charged fermions, the new 
CP phase has to be about 0.1 or less. An ongoing electron EDM experiment will improve the electron 
EDM bound by about one order of magnitude in the near future, potentially worsening the Higgs CP 
problem. In models with charged scalars to enhance Higgs diphoton rate, the two-loop EDM could 
be negligible if there is no light pseudoscalar mixing with the Higgs. Thus, in the case of charged 
scalars, the Higgs CP problem would be less robust. However, vacuum instability problems constrain 
the charged scalar models severely. 

It is remarkable that experiments underground (direct detection), in atomic physics laboratories 
(EDM measurements), in colliders, and in space (Fermi- LAT) are setting comparable bounds on new 
charged matter. This offers the prospect of a stunning coincidence of signals, if such matter exists. 
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A Fermionic DM through a box diagram 

In this appendix we will summarize some results on one further model, to show that the qualitative 
conclusions reached in Section 2.1 are not strongly dependent on the initial assumption that all particles 
involved are scalars. We assume DM is a fermion £, £ c annihilating through loops of charged fermions 
ip,ip c an d a charged scalar <fi: 

m l H 2 + (3?^ c f + fff^r^ + ™?f ? + m xX c X + m^ c tp + yH^ X c + J/ c #Vx + h.c.) . (A.l) 
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The fermions ift, \ are, as before, a singlet and doublet, while (ft is a charge 1 scalar. In this case the 
RGEs are: 



16tH /%) = 2 5 f + gf g ( +y CJ g ( , 

i6n 2 p(y) = ~y 3 + v c2 y + gfy, 



(A.2) 
(A.3) 



together with the analogous equations where couplings with and without superscript c are inter- 
changed. We fix the masses my, m x and the Yukawas y,y c as in Section 2.2 to give a 50% en- 
hancement of the h —¥ 77 rate. In this case we don't need to include a CP-violating phase, so we 
take m x — my — 335 GeV, y = y c = 1.12, which produces charged fermion mass eigenstates at 
fn\ t i = 140, 530 GeV. The mass matrix for the fermions is diagonalized by two matrices U and V, i.e. 
U*A4V> is the diagonal matrix with eigenvalues m\ and m^ (as in, e.g., Ref. [86]). 

We fix = 130 GeV. For simplicity we will fix g| 
loop of charged fermions and scalars is given by [25, 79, 



gl. Then the annihilation rate through a 



2 2 

a m| 



647r 3 r 



E2„ T , "> 



to 2 +m 2 





2TO,?Tii 



(A.4) 



with the loop functions defined in Ref. [79]. Then we can achieve a dark matter annihilation to 
diphoton rate of 10 -27 cm 3 /s with the choice: = 150 GeV, — 3.6. As in the case of an all-scalar 
loop, this is a rather large coupling; in fact, with our definition of normalized coupling from Sec. 2.1.4, 
solving the RGE shows that g^ and g| reach a normalized value of 1 already at a renormalization scale 
of 235 GeV. They reach a value of 4ir well below a TeV. Hence, in the case of a box diagram, one 
really should imagine that the theory is one of composite dark matter coupling to composite charged 
fields nearby in mass, much as in the all-scalar case we considered above. 

In this case there is again a loop- induced dark matter coupling to the Higgs. We have not fully 
calculated the diagram, but schematically one expects a coupling of the form 



g^yy m x 



16tt 2 A 2 



(A.5) 



with A of order the masses of particles running in the loop. Taking A ~ 150 GeV and plugging in our 
choices, we find a coupling ~ 0.4/i££ c , in conflict with direct detection bounds unless the numerical 
coefficient proves to be quite small. We expect that a more careful calculation would continue to 
support the same qualitative story as in Section 2.1: without resonant enhancement, the couplings 
involved in fitting the line rate and a large h — > 77 excess are so large as to pose a strong tension with 
the absence of large direct dark matter couplings to the Higgs implied by XENON. 



B Counting CP-odd operators 

The literature can be slightly confusing on the enumeration of the set of dimension-six CP-odd op- 
erators constructed from electroweak gauge bosons and the Higgs. In particular, one can find two 
more CP-odd dimension six operators besides the ones we listed in Sec. 3.1: (D ^H)^ ' o a D v HW^ v and 
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[PfiHy D V HB^ [91, 92]. We will show that these two operators are related by equations of motion 
to the ones we listed. We will take {D^Hy a a D u HW^ v for example. 

{D ll H)^a a D v UW^ a = (d^ + i 9 -H^W b y + i 9 -H^ B^Y)a a D v HW^ a 
= -H^a a D v Hd ll W^ a - iJ t a a (d^D u H)W^ a + (i^H^W b a b + i^H^B^y) a a D u HW fn '' a 
= -H^o- a D iy Hd fl W^ a _ H ] a a Uq^ + il W b a b + i^B^ D„H^J W^ a + gHh abc W b <r c D„HW^' a 
= -H^a a D u H(-ge abc W b W^ c ) - H^a a (D^D u H)W^ a + gHU abc W b <r c D„HW^ a 
= -^a a [D^,D u ]HW^ a 

= ~Hia a (gW*y + g'B^HW^ = '-gH^HW^W^ + -^g' o a HB^W^ a , (B.l) 

where in the second line we integrated by parts; in the third line we used [<7 a ,cr fc ] = 2ie abc a c ; in the 
third line we used Bianchi identity e^ pa (D tl W pa ) a = 0; in the last line we used a a a b = 5 ab I + ie abc a c . 
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